Functional proteomics of colon cancer Consensus Molecular Subtypes.

BACKGROUND: The Colorectal Cancer Subtyping Consortium established four Consensus Molecular Subtypes (CMS) in colorectal cancer: CMS1 (microsatellite-instability [MSI], Immune), CMS2 (Canonical, epithelial), CMS3 (Metabolic), and CMS4 (Mesenchymal). However, only MSI tumour patients have seen a change in their disease management in clinical practice. This study aims to characterise the proteome of colon cancer CMS and broaden CMS's clinical utility. METHODS: One-hundred fifty-eight paraffin samples from stage II-III colon cancer patients treated with adjuvant chemotherapy were analysed through DIA-based mass-spectrometry proteomics. RESULTS: CMS1 exhibited overexpression of immune-related proteins, specifically related to neutrophils, phagocytosis, antimicrobial response, and a glycolytic profile. These findings suggested potential therapeutic strategies involving immunotherapy and glycolytic inhibitors. CMS3 showed overexpression of metabolic proteins. CMS2 displayed a heterogeneous protein profile. Notably, two proteomics subtypes within CMS2, with different protein characteristics and prognoses, were identified. CMS4 emerged as the most distinct group, featuring overexpression of proteins related to angiogenesis, extracellular matrix, focal adhesion, and complement activation. CMS4 showed a high metastatic profile and suggested possible chemoresistance that may explain its worse prognosis. CONCLUSIONS: DIA proteomics revealed new features for each colon cancer CMS subtype. These findings provide valuable insights into potential therapeutic targets for colorectal cancer subtypes in the future.

Multi-omics Characterization of Response to PD-1 Inhibitors in Advanced Melanoma.

Immunotherapy improves the survival of patients with advanced melanoma, 40% of whom become long-term responders. However, not all patients respond to immunotherapy. Further knowledge of the processes involved in the response and resistance to immunotherapy is still needed. In this study, clinical paraffin samples from fifty-two advanced melanoma patients treated with anti-PD-1 inhibitors were assessed via high-throughput proteomics and RNA-seq. The obtained proteomics and transcriptomics data were analyzed using multi-omics network analyses based on probabilistic graphical models to identify those biological processes involved in the response to immunotherapy. Additionally, proteins related to overall survival were studied. The activity of the node formed by the proteins involved in protein processing in the endoplasmic reticulum and antigen presentation machinery was higher in responders compared to non-responders; the activity of the immune and inflammatory response node was also higher in those with complete or partial responses. A predictor for overall survival based on two proteins (AMBP and PDSM5) was defined. In summary, the response to anti-PD-1 therapy in advanced melanoma is related to protein processing in the endoplasmic reticulum, and also to genes involved in the immune and inflammatory responses. Finally, a two-protein predictor can define survival in advanced disease. The molecular characterization of the mechanisms involved in the response and resistance to immunotherapy in melanoma leads the way to establishing therapeutic alternatives for patients who will not respond to this treatment.

Metabolomics unravels subtype-specific characteristics related to neoadjuvant therapy response in breast cancer patients.

INTRODUCTION: Breast cancer is the most diagnosed tumor and the leading cause of cancer death in women worldwide. Metabolomics allows the quantification of the entire set of metabolites in blood samples, making it possible to study differential metabolomics patterns related to neoadjuvant treatment in the breast cancer neoadjuvant setting. OBJECTIVES: Characterizing metabolic differences in breast cancer blood samples according to their response to neoadjuvant treatment. METHODS: One hundred and three plasma samples of breast cancer patients, before receiving neoadjuvant treatment, were analyzed through UPLC-MS/MS metabolomics. Then, metabolomics data were analyzed using probabilistic graphical models and biostatistics methods. RESULTS: Metabolomics data allowed the identification of differences between groups according to response to neoadjuvant treatment. These differences were specific to each breast cancer subtype. Patients with HER2+ tumors showed differences in metabolites related to amino acids and carbohydrates pathways between the two pathological response groups. However, patients with triple-negative tumors showed differences in metabolites related to the long-chain fatty acids pathway. Patients with Luminal B tumors showed differences in metabolites related to acylcarnitine pathways. CONCLUSIONS: It is possible to identify differential metabolomics patterns between complete and partial responses to neoadjuvant therapy, being this metabolomic profile specific for each breast cancer subtype.

Utility of CYP2D6 copy number variants as prognostic biomarker in localized anal squamous cell carcinoma.

BACKGROUND: Anal squamous cell carcinoma (ASCC) is an infrequent tumor whose treatment has not changed since the 1970s. The aim of this study is the identification of biomarkers allowing personalized treatments and improvement of therapeutic outcomes. METHODS: Forty-six paraffin tumor samples from ASCC patients were analyzed by whole-exome sequencing. Copy number variants (CNVs) were identified and their relation to disease-free survival (DFS) was studied and validated in an independent retrospective cohort of 101 ASCC patients from the Multidisciplinary Spanish Digestive Cancer Group (GEMCAD). GEMCAD cohort proteomics allowed assessing the biological features of these tumors. RESULTS: On the discovery cohort, the median age was 61 years old, 50% were males, stages I/II/III: 3 (7%)/16 (35%)/27 (58%), respectively, median DFS was 33 months, and overall survival was 45 months. Twenty-nine genes whose duplication was related to DFS were identified. The most representative was duplications of the CYP2D locus, including CYP2D6, CYP2D7P, and CYP2D8P genes. Patients with CYP2D6 CNV had worse DFS at 5 years than those with two CYP2D6 copies (21% vs. 84%; p < .0002, hazard ratio [HR], 5.8; 95% confidence interval [CI], 2.7-24.9). In the GEMCAD validation cohort, patients with CYP2D6 CNV also had worse DFS at 5 years (56% vs. 87%; p = .02, HR = 3.6; 95% CI, 1.1-5.7). Mitochondria and mitochondrial cell-cycle proteins were overexpressed in patients with CYP2D6 CNV. CONCLUSIONS: Tumor CYP2D6 CNV identified patients with a significantly worse DFS at 5 years among localized ASCC patients treated with 5-fluorouracil, mitomycin C, and radiotherapy. Proteomics pointed out mitochondria and mitochondrial cell-cycle genes as possible therapeutic targets for these high-risk patients. PLAIN LANGUAGE SUMMARY: Anal squamous cell carcinoma is an infrequent tumor whose treatment has not been changed since the 1970s. However, disease-free survival in late staged tumors is between 40% and 70%. The presence of an alteration in the number of copies of CYP2D6 gene is a biomarker of worse disease-free survival. The analysis of the proteins in these high-risk patients pointed out mitochondria and mitochondrial cell-cycle genes as possible therapeutic targets. Therefore, the determination of the number of copies of CYP2D6 allows the identification of anal squamous carcinoma patients with a high-risk of relapse that could be redirected to a clinical trial. Additionally, this study may be useful to suggest new treatment strategies to increase current therapy efficacy.

Sorting Transcriptomics Immune Information from Tumor Molecular Features Allows Prediction of Response to Anti-PD1 Therapy in Patients with Advanced Melanoma.

Immunotherapy based on anti-PD1 antibodies has improved the outcome of advanced melanoma. However, prediction of response to immunotherapy remains an unmet need in the field. Tumor PD-L1 expression, mutational burden, gene profiles and microbiome profiles have been proposed as potential markers but are not used in clinical practice. Probabilistic graphical models and classificatory algorithms were used to classify melanoma tumor samples from a TCGA cohort. A cohort of patients with advanced melanoma treated with PD-1 inhibitors was also analyzed. We established that gene expression data can be grouped in two different layers of information: immune and molecular. In the TCGA, the molecular classification provided information on processes such as epidermis development and keratinization, melanogenesis, and extracellular space and membrane. The immune layer classification was able to distinguish between responders and non-responders to immunotherapy in an independent series of patients with advanced melanoma treated with PD-1 inhibitors. We established that the immune information is independent than molecular features of the tumors in melanoma TCGA cohort, and an immune classification of these tumors was established. This immune classification was capable to determine what patients are going to respond to immunotherapy in a new cohort of patients with advanced melanoma treated with PD-1 inhibitors Therefore, this immune signature could be useful to the clinicians to identify those patients who will respond to immunotherapy.

Proteomics Characterization of Clear Cell Renal Cell Carcinoma.

PURPOSE: To explore the tumor proteome of patients diagnosed with localized clear cell renal cancer (ccRCC) and treated with surgery. MATERIAL AND METHODS: A total of 165 FFPE tumor samples from patients diagnosed with ccRCC were analyzed using DIA-proteomics. Proteomics ccRCC subtypes were defined using a consensus cluster algorithm (CCA) and characterized by a functional approach using probabilistic graphical models and survival analyses. RESULTS: We identified and quantified 3091 proteins, including 2026 high-confidence proteins. Two proteomics subtypes of ccRCC (CC1 and CC2) were identified by CC using the high-confidence proteins only. Characterization of molecular differences between CC1 and CC2 was performed in two steps. First, we defined 514 proteins showing differential expression between the two subtypes using a significance analysis of microarrays analysis. Proteins overexpressed in CC1 were mainly related to translation and ribosome, while proteins overexpressed in CC2 were mainly related to focal adhesion and membrane. Second, a functional analysis using probabilistic graphical models was performed. CC1 subtype is characterized by an increased expression of proteins related to glycolysis, mitochondria, translation, adhesion proteins related to cytoskeleton and actin, nucleosome, and spliceosome, while CC2 subtype showed higher expression of proteins involved in focal adhesion, extracellular matrix, and collagen organization. CONCLUSIONS: ccRCC tumors can be classified in two different proteomics subtypes. CC1 and CC2 present specific proteomics profiles, reflecting alterations of different molecular pathways in each subtype. The knowledge generated in this type of studies could help in the development of new drugs targeting subtype-specific deregulated pathways.

MiRNAs in renal cell carcinoma

MicroRNAs (miRNAs) are small RNA sequences that act as post-transcriptional regulatory genes to control many cellular processes through pairing bases with a complementary messenger RNA (mRNA). A single miRNA molecule can regulate more than 200 different transcripts and the same mRNA can be regulated by multiple miRNAs. In this review, we highlight the importance of miRNAs and collect the existing evidence on their relationship with kidney cancer. (AU)

MiRNAs in renal cell carcinoma.

MicroRNAs (miRNAs) are small RNA sequences that act as post-transcriptional regulatory genes to control many cellular processes through pairing bases with a complementary messenger RNA (mRNA). A single miRNA molecule can regulate more than 200 different transcripts and the same mRNA can be regulated by multiple miRNAs. In this review, we highlight the importance of miRNAs and collect the existing evidence on their relationship with kidney cancer.

Identification of Carcinogenesis and Tumor Progression Processes in Pancreatic Ductal Adenocarcinoma Using High-Throughput Proteomics.

Pancreatic ductal adenocarcinoma (PDAC) is an aggressive disease with an overall 5-year survival rate of just 5%. A better understanding of the carcinogenesis processes and the mechanisms of the progression of PDAC is mandatory. Fifty-two PDAC patients treated with surgery and adjuvant therapy, with available primary tumors, normal tissue, preneoplastic lesions (PanIN), and/or lymph node metastases, were selected for the study. Proteins were extracted from small punches and analyzed by LC-MS/MS using data-independent acquisition. Proteomics data were analyzed using probabilistic graphical models, allowing functional characterization. Comparisons between groups were made using linear mixed models. Three proteomic tumor subtypes were defined. T1 (32% of patients) was related to adhesion, T2 (34%) had metabolic features, and T3 (34%) presented high splicing and nucleoplasm activity. These proteomics subtypes were validated in the PDAC TCGA cohort. Relevant biological processes related to carcinogenesis and tumor progression were studied in each subtype. Carcinogenesis in the T1 subtype seems to be related to an increase of adhesion and complement activation node activity, whereas tumor progression seems to be related to nucleoplasm and translation nodes. Regarding the T2 subtype, it seems that metabolism and, especially, mitochondria act as the motor of cancer development. T3 analyses point out that nucleoplasm, mitochondria and metabolism, and extracellular matrix nodes could be involved in T3 tumor carcinogenesis. The identified processes were different among proteomics subtypes, suggesting that the molecular motor of the disease is different in each subtype. These differences can have implications for the development of future tailored therapeutic approaches for each PDAC proteomics subtype.

Molecular characterization of triple negative breast cancer formaldehyde-fixed paraffin-embedded samples by data-independent acquisition proteomics.

Triple negative breast cancer accounts for 15%-20% of all breast carcinomas and is clinically characterized by an aggressive phenotype and poor prognosis. Triple negative tumors do not benefit from targeted therapies, so further characterization is needed to define subgroups with potential therapeutic value. In this work, the proteomes of 125 formalin-fixed paraffin-embedded samples from patients diagnosed with non-metastatic triple negative breast cancer were analyzed using data-independent acquisition + in a LTQ-Orbitrap Fusion Lumos mass spectrometer coupled to an EASY-nLC 1000. 1206 proteins were identified in at least 66% of the samples. Hierarchical clustering, probabilistic graphical models and Significance Analysis of Microarrays were combined to characterize proteomics-based molecular groups. Two molecular groups were defined with differences in biological processes such as glycolysis, translation and immune response. These two molecular groups showed also several differentially expressed proteins. This clinically homogenous dataset may serve to design new therapeutic strategies in the future.

Protein content of blood-derived extracellular vesicles: An approach to the pathophysiology of cerebral hemorrhage.

Introduction: Extracellular vesicles (EVs) participate in cell-to-cell paracrine signaling and can be biomarkers of the pathophysiological processes underlying disease. In intracerebral hemorrhage, the study of the number and molecular content of circulating EVs may help elucidate the biological mechanisms involved in damage and repair, contributing valuable information to the identification of new therapeutic targets. Methods: The objective of this study was to describe the number and protein content of blood-derived EVs following an intracerebral hemorrhage (ICH). For this purpose, an experimental ICH was induced in the striatum of Sprague-Dawley rats and EVs were isolated and characterized from blood at baseline, 24 h and 28 days. The protein content in the EVs was analyzed by mass spectrometric data-dependent acquisition; protein quantification was obtained by sequential window acquisition of all theoretical mass spectra data and compared at pre-defined time points. Results: Although no differences were found in the number of EVs, the proteomic study revealed that proteins related to the response to cellular damage such as deubiquitination, regulation of MAP kinase activity (UCHL1) and signal transduction (NDGR3), were up-expressed at 24 h compared to baseline; and that at 28 days, the protein expression profile was characterized by a higher content of the proteins involved in healing and repair processes such as cytoskeleton organization and response to growth factors (COR1B) and the regulation of autophagy (PI42B). Discussion: The protein content of circulating EVs at different time points following an ICH may reflect evolutionary changes in the pathophysiology of the disease.

Description of the genetic variants identified in a cohort of patients diagnosed with localized anal squamous cell carcinoma and treated with panitumumab.

Squamous cell carcinoma is the most frequent histologic type of anal carcinoma. The standard of care since the 1970s has been a combination of 5-fluorouracil, mitomycin C, and radiotherapy. This treatment is very effective in T1/T2 tumors (achieving complete regression in 80-90% of tumors). However, in T3/T4 tumors, the 3-year relapse free survival rate is only 50%. The VITAL trial aimed to assess the efficacy and safety of panitumumab in combination with this standard treatment. In this study, 27 paraffin-embedded samples from the VITAL trial and 18 samples from patients from daily clinical practice were analyzed by whole-exome sequencing and the influence of the presence of genetic variants in the response to panitumumab was studied. Having a moderate- or high-impact genetic variant in PIK3CA seemed to be related to the response to panitumumab. Furthermore, copy number variants in FGFR3, GRB2 and JAK1 were also related to the response to panitumumab. These genetic alterations have also been studied in the cohort of patients from daily clinical practice (not treated with panitumumab) and they did not have a predictive value. Therefore, in this study, a collection of genetic alterations related to the response with panitumumab was described. These results could be useful for patient stratification in new anti-EGFR clinical trials.

Bayesian networks established functional differences between breast cancer subtypes.

Breast cancer is a heterogeneous disease. In clinical practice, tumors are classified as hormonal receptor positive, Her2 positive and triple negative tumors. In previous works, our group defined a new hormonal receptor positive subgroup, the TN-like subtype, which had a prognosis and a molecular profile more similar to triple negative tumors. In this study, proteomics and Bayesian networks were used to characterize protein relationships in 96 breast tumor samples. Components obtained by these methods had a clear functional structure. The analysis of these components suggested differences in processes such as mitochondrial function or extracellular matrix between breast cancer subtypes, including our new defined subtype TN-like. In addition, one of the components, mainly related with extracellular matrix processes, had prognostic value in this cohort. Functional approaches allow to build hypotheses about regulatory mechanisms and to establish new relationships among proteins in the breast cancer context.

Computational models applied to metabolomics data hints at the relevance of glutamine metabolism in breast cancer.

BACKGROUND: Metabolomics has a great potential in the development of new biomarkers in cancer and it has experiment recent technical advances. METHODS: In this study, metabolomics and gene expression data from 67 localized (stage I to IIIB) breast cancer tumor samples were analyzed, using (1) probabilistic graphical models to define associations using quantitative data without other a priori information; and (2) Flux Balance Analysis and flux activities to characterize differences in metabolic pathways. RESULTS: On the one hand, both analyses highlighted the importance of glutamine in breast cancer. Moreover, cell experiments showed that treating breast cancer cells with drugs targeting glutamine metabolism significantly affects cell viability. On the other hand, these computational methods suggested some hypotheses and have demonstrated their utility in the analysis of metabolomics data and in associating metabolomics with patient's clinical outcome. CONCLUSIONS: Computational analyses applied to metabolomics data suggested that glutamine metabolism is a relevant process in breast cancer. Cell experiments confirmed this hypothesis. In addition, these computational analyses allow associating metabolomics data with patient prognosis.

miRNA profiling in renal carcinoma suggest the existence of a group of pro-angionenic tumors in localized clear cell renal carcinoma.

Renal cell carcinoma comprises a variety of entities, the most common being the clear-cell, papillary and chromophobe subtypes. These subtypes are related to different clinical evolution; however, most therapies have been developed for clear-cell carcinoma and there is not a specific treatment based on different subtypes. In this study, one hundred and sixty-four paraffin samples from primary nephrectomies for localized tumors were analyzed. MiRNAs were isolated and measured by microRNA arrays. Significance Analysis of Microarrays and Consensus Cluster algorithm were used to characterize different renal subtypes. The analyses showed that chromophobe renal tumors are a homogeneous group characterized by an overexpression of miR 1229, miR 10a, miR 182, miR 1208, miR 222, miR 221, miR 891b, miR 629-5p and miR 221-5p. On the other hand, clear cell renal carcinomas presented two different groups inside this histological subtype, with differences in miRNAs that regulate focal adhesion, transcription, apoptosis and angiogenesis processes. Specifically, one of the defined groups had an overexpression of proangiogenic microRNAs miR185, miR126 and miR130a. In conclusion, differences in miRNA expression profiles between histological renal subtypes were established. In addition, clear cell renal carcinomas had different expression of proangiogenic miRNAs. With the emergence of antiangiogenic drugs, these differences could be used as therapeutic targets in the future or as a selection method for tailoring personalized treatments.

Genetic Profile and Functional Proteomics of Anal Squamous Cell Carcinoma: Proposal for a Molecular Classification.

Anal squamous cell carcinoma is a rare tumor. Chemo-radiotherapy yields a 50% 3-year relapse-free survival rate in advanced anal cancer, so improved predictive markers and therapeutic options are needed. High-throughput proteomics and whole-exome sequencing were performed in 46 paraffin samples from anal squamous cell carcinoma patients. Hierarchical clustering was used to establish groups de novo Then, probabilistic graphical models were used to study the differences between groups of patients at the biological process level. A molecular classification into two groups of patients was established, one group with increased expression of proteins related to adhesion, T lymphocytes and glycolysis; and the other group with increased expression of proteins related to translation and ribosomes. The functional analysis by the probabilistic graphical model showed that these two groups presented differences in metabolism, mitochondria, translation, splicing and adhesion processes. Additionally, these groups showed different frequencies of genetic variants in some genes, such as ATM, SLFN11 and DST Finally, genetic and proteomic characteristics of these groups suggested the use of some possible targeted therapies, such as PARP inhibitors or immunotherapy.

Computational metabolism modeling predicts risk of distant relapse-free survival in breast cancer patients.

Aim: Differences in metabolism among breast cancer subtypes suggest that metabolism plays an important role in this disease. Flux balance analysis is used to explore these differences as well as drug response. Materials & methods: Proteomics data from breast tumors were obtained by mass-spectrometry. Flux balance analysis was performed to study metabolic networks. Flux activities from metabolic pathways were calculated and used to build prognostic models. Results: Flux activities of vitamin A, tetrahydrobiopterin and ß-alanine metabolism pathways split our population into low- and high-risk patients. Additionally, flux activities of glycolysis and glutamate metabolism split triple negative tumors into low- and high-risk groups. Conclusion: Flux activities summarize flux balance analysis data and can be associated with prognosis in cancer.

Melanoma proteomics suggests functional differences related to mutational status.

Melanoma is the most lethal cutaneous cancer. New drugs have recently appeared; however, not all patients obtain a benefit of these new drugs. For this reason, it is still necessary to characterize melanoma at molecular level. The aim of this study was to explore the molecular differences between melanoma tumor subtypes, based on BRAF and NRAS mutational status. Fourteen formalin-fixed, paraffin-embedded melanoma samples were analyzed using a high-throughput proteomics approach, combined with probabilistic graphical models and Flux Balance Analysis, to characterize these differences. Proteomics analyses showed differences in expression of proteins related with fatty acid metabolism, melanogenesis and extracellular space between BRAF mutated and BRAF non-mutated melanoma tumors. Additionally, probabilistic graphical models showed differences between melanoma subgroups at biological processes such as melanogenesis or metabolism. On the other hand, Flux Balance Analysis predicts a higher tumor growth rate in BRAF mutated melanoma samples. In conclusion, differential biological processes between melanomas showing a specific mutational status can be detected using combined proteomics and computational approaches.

Isotopologue Multipoint Calibration for Proteomics Biomarker Quantification in Clinical Practice.

Targeted proteomics has become the method of choice for biomarker validation in human biopsies due to its high sensitivity, reproducibility, accuracy, and precision. However, for targeted proteomics to be transferred to clinical routine there is the need to reduce its complexity, make its procedures simpler, increase its throughput, and improve its analytical performance. Here we present the Isotopologue Multipoint Calibration (ImCal) quantification strategy, which uses a mix of isotopologue peptides to generate internal multipoint calibration curves for each individual sample and to accurately quantify biomarker peptides in clinical applications without the need of expert supervision. ImCal relies on the use of five different isotopically-labelled peptides of different nominal mass mixed at different concentrations to be used as an internal calibration curve for each endogenous peptide. The use of internal multipoint calibration curves is well-suited for the generation of ready-to-use biomarker kits for clinical applications as it is compatible with both high- and low-resolution mass spectrometers and different levels of endogenous peptide, it eliminates the need for blank matrixes required in external curves, it allows the evaluation of matrix effects and the valid quantification range in each individual sample, and it does not require expert adjustment. We used the ImCal method to quantify HER2 in 35 breast cancer formalin-fixed paraffin-embedded patient samples, revealing a high degree of heterogeneity among patients, which contrasts with the homogeneous immunohistochemistry patient classification. Our work illustrates how an improvement of mass spectrometry methods for biomarker quantification can provide fine-grain patient stratification, and thus better disease diagnostic and prognosis.

A novel approach to triple-negative breast cancer molecular classification reveals a luminal immune-positive subgroup with good prognoses.

Triple-negative breast cancer is a heterogeneous disease characterized by a lack of hormonal receptors and HER2 overexpression. It is the only breast cancer subgroup that does not benefit from targeted therapies, and its prognosis is poor. Several studies have developed specific molecular classifications for triple-negative breast cancer. However, these molecular subtypes have had little impact in the clinical setting. Gene expression data and clinical information from 494 triple-negative breast tumors were obtained from public databases. First, a probabilistic graphical model approach to associate gene expression profiles was performed. Then, sparse k-means was used to establish a new molecular classification. Results were then verified in a second database including 153 triple-negative breast tumors treated with neoadjuvant chemotherapy. Clinical and gene expression data from 494 triple-negative breast tumors were analyzed. Tumors in the dataset were divided into four subgroups (luminal-androgen receptor expressing, basal, claudin-low and claudin-high), using the cancer stem cell hypothesis as reference. These four subgroups were defined and characterized through hierarchical clustering and probabilistic graphical models and compared with previously defined classifications. In addition, two subgroups related to immune activity were defined. This immune activity showed prognostic value in the whole cohort and in the luminal subgroup. The claudin-high subgroup showed poor response to neoadjuvant chemotherapy. Through a novel analytical approach we proved that there are at least two independent sources of biological information: cellular and immune. Thus, we developed two different and overlapping triple-negative breast cancer classifications and showed that the luminal immune-positive subgroup had better prognoses than the luminal immune-negative. Finally, this work paves the way for using the defined classifications as predictive features in the neoadjuvant scenario.